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Abstract Accurate fitting formulae to the synchrotron function, F{x), and its complementary 

■ function, G{x), are performed and presented. The corresponding relative errors are less than 
OO ■ 

' 0.26% and 0.035% for F{x) and G{x), respectively. To this aim we have, first, fitted the 

G^ ! 

\^ , modified Bessel functions, K^/^{x) and K2/j,{x). For all the fitted functions, the general fit 

expression is the same, and is based on the well known asymptotic forms for low and large 

o 

■ x-values for each function. It consists of multiplying each asymptotic form by a function 
that tends to unity or zero for low and large a--values. Simple formulae are suggested in 
this paper, depending on adjustable parameters. The latter have been determined by adopting 

. the Levenberg-Marquardt algorithm. The proposed formulae should be of great utility and 



simplicity for computing spectral powers and the degree of polarization for the synchrotron 
radiation, both for laboratory and astrophysical appUcations. 
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1 INTRODUCTION 

Analytical approximate formulae are often very useful and may be indispensable in order to avoid the 
computation of complicated transcendental functions. This is the case of the modified Bessel functions and 
their integrals, especially those of the second kind with fractional order, e.g., K<:^/^{x) and K2/z{x), on 
which we focus our attention in this contribution. We start by presenting, in section 2, results of fits to these 
two functions. Then, in section , we deduce the expression of the complementary synchrotron function, 
G{x) = xK2/3,{x), directly from function K2/3{x), and report the corresponding fit to the synchrotron 
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2 MODIFIED BESSEL FUNCTIONS AND K2/3 
2.1 Definitions 

The modified Bessel functions, I±„{x) and K^{x), of the first and second kind, respectively, are particular 
solutions of Bessel's cylindrical differential equation, i.e., (Abramowitz & Stegun 1965) 

x^--^+x- {x^ + v^)w^Q. (1) 

aa;^ ax ^ ' 

Function K^{x) expresses as (Abramowitz & Stegun 1965) 

, , TT I_,,(x) — I, Ax) 

Ku{x) = tt^^^tV^^, (2) 
2 sm(i^7r) 

in terms of function Ii,{x) that writes as (Abramowitz & Stegun 1965) 

in form of an ascending series involving the F function. Besides, function K„{x) can also be written as 
(Abramowitz & Stegun 1965) 

K,{x) = / e-^-\e - ly-Ut, (4) 



in integral representation. 

Finally, this function admits the following simplified asymptotic forms (Abramowitz & Stegun 1965) 



A.(..)4r(.)(|) 



for a; <C 1 



A2{x) = J—x 2e ^ fora;>l 



2.2 Fitting formulae 



In fitting a function, f{x) (here, the modified Bessel functions and the synchrotron functions), the main idea 
consists in expressing it in terms of its known asymptotic forms, say Ai (x) for low x-values and A2 (x) for 
large z-values, and to put it under the form 

fix)=Ai{x)6iix)+A2{x)S2{x), (6) 
where 5i (x) and 82 (x) are the functions one is looking for, which must respectively obey the limits 

61 (x) w 1 for a; < 1 

(7) 

Si (x) w for X- > 1 

and 



S2{x) « for a; < 1 

(8) 
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Table 1 Coefficients a^^' and a^^-* for function 1^5/3(2;). 



k 


4.^' 


4^' 


1 


-1.0194198041210243 


-15.761577796582387 


2 


+0.28011396300530672 




3 


-7.71058491739234908 x 10"^ 





Notes: With this set of coefficients, the relative error is < 0.48%. 

Table 2 Coefficients a^^-* and a^^'' for function K2/3{x). 



k 




(2) 


1 


-1.3746667760953621 


-0.33550751062084 


2 


+0.44040512552162292 




3 


-0.15527012012316799 





Notes; With this set of coefficients, the relative error is < 0.54%. 

Table 3 Coefficients a^^-* and a^^^ for function K2/3{x). 

1 -1.0010216415582440 -0.2493940736333195 

2 +0.88350305221249859 +0.9122693061687756 

3 -3.6240174463901829 +1.2051408667145216 

4 +0.57393980442916881 -5.5227048291651126 

Notes: With this set of coefficients, the relative error is < 0.035%. 



For this purpose, we propose the following expressions: 



Til 



(1)3,1//= 



fc=l 



and 



S2{X) l_e^2W 



(2)2,l/fc 



fe=l 



(9) 



(10) 



In order to extract coefficients a), ' and aj, for a given couple of orders (711,712), we proceed by 
chi-squares minimization with adopting the Levenberg-Marquardt algorithm (Levenberg 1944, Marquardt 
1963), in log-log scale. The obtained fit results to functions Kr,/^{x) and K2/3{x) are presented in tables 1 
and 2, respectively, in terms of coefficients a^^^ and a^^\ with 711 = 3 and 712 = 1 and relative respective 
errors, < 0.48% and < 0.54%. These fits to functions 1^5/3(2;) and -^^2/3(2^) are plotted in figures 1 and 3, 
respectively, while the corresponding relative errors are reported in figures 2 and 4. 

For high accuracy, we give, in table 3, fit results for function K2/3{x), with ni = 712 = 4 and with a 
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X 

Fig. 1 The modified Bessel function, K^/^{x), together with its corresponding fit according to 
equation (6). 
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Fig. 2 The relative error for the modified Bessel function, K^/^^lx), corresponding to the set of 
coefficients reported by table 1 . 



3 SYNCHROTRON FUNCTIONS 
3.1 Definitions 

The synchrotron functions, F{x) and G(x), are defined by ("Westfold 1959; Jackson 1962; Rybicki & 
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Fig. 3 The modified Bessel function, K2/^{x), together with its corresponding fit, according to 
equation (6). 



Relative Error for Modified Bessel function K ^(x) 
0.006 ' . 
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Fig.4 The relative error for the modified Bessel function, K2/^{x), corresponding to the set of 
coefficients reported by table 2. 



F{x) = X / K^,^{x')dx' 

J'- . (11) 

G{x) = xK2/3ix) 

Function G{x) is called the complementary synchrotron function and is sometimes noted Fp{x) 
(Westfold 1959). The corresponding simplest asymptotic forms of these functions have the following ex- 
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Table 4 Coefficients a^^"* and a^^^ for the synchrotron function 

F{x). 



k 


4^' 




1 


-0.97947838884478688 


-4.69247165562628882 x 10"^ 


2 


-0.83333239129525072 


-0.70055018056462881 


3 


+0.15541796026816246 


1.03876297841949544 x 10"^ 



Notes: With this set of coefficients, the relative error is < 0.26%. 




Fig. 5 The synchrotron function, F{x), together with its corresponding fit according to equation 
(6). 



Fix) 



Fix for a; < 1 
F2e~''x^/^ for x > 1 



(12) 



and 



Gix) 



Gix^^^ for a: < 1 
Gae-'^.T^/^ for X > 1 ' 

where Fi = 7r25/3/\/3r(l/3), F2 = and d = Fi/2 , G2 = F2. 



(13) 



3.2 Fitting formulae 

Function G{x) can be easily derived directly from the fit to function K2/3{x). One has just to multiply 
the latter by variable x. For fitting function F{x) , we proceed in the same way as for the modified Bessel 
functions, i.e., putting it under the form given by equation (6). We have just to consider the corresponding 
asymptotic forms given by equation (12). The corresponding fit coefficients are reported in table 4. With 
these coefficients, the relative error is < 0.26%. Function F{x) is plotted in figure 5, together with the 



Analytical Fits to tlie Synclirotron Functions 



0.004 



Relative Error for Synclirotron function F{x) 
"T — I — I — I — I — I — I 



8 

w 



0.002 



0.000 



-0.002 



-0.004 I I 

to' 10^ lO' lO' lO' 10° to' lO' 

X 

Fig. 6 The relative error for the synchrotron function, F{x), corresponding to the set of coeffi- 
cients reported by table 4. 

4 CONCLUSION 

We have presented analytical fit formulae with good accuracies for the synchrotron function, F{x), and 
its complementary function, G{x), based on their known asymptotic forms for low and large x-values. We 
propose these formulae for the aim of directly and simply computing these transcendental functions with 
avoiding fastidious calculations. The derived general fit formulae can thus be used to evaluate the modified 
Bessel functions of any order: integer or non integer Finally, these fit formulae should be of great help for 
computing quantities of interest to synchrotron radiation such as, e.g., the spectral power and the degree of 
polarization, both for laboratory and astrophysical applications. 
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